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ABSTRACT 

We present the results from simultaneous observations of the Narrow-Line Seyfert 1 galaxy NGC 405 1 with the 
Chandra High Energy Transmission Grating Spectrometer and the HST Space Telescope Imaging Spectrograph. 
The X-ray grating spectrum reveals absorption and emission lines from hydrogen-like and helium-like ions of O, 
Ne, Mg and Si. We resolve two distinct X-ray absorption systems: a high-velocity blueshifted system at -2340 ± 
130 km s" 1 and a low-velocity blueshifted system at -600 ±130 km s . In the UV spectrum we detect strong 
absorption, mainly from C IV, N V and Si IV, that is resolved into as many as nine different intrinsic absorption 
systems with velocities between -650 km s" 1 and 30 km s" 1 . Although the low-velocity X-ray absorption is 
consistent in velocity with many of the UV absorption systems, the high-velocity X-ray absorption seems to have 
no UV counterpart. In addition to the absorption and emission lines, we also observe rapid X-ray variability and a 
state of low X-ray flux during the last w 15 ks of the observation. NGC 405 1 has a soft X-ray excess which we fit 
in both the high and low X-ray flux states. The high-resolution X-ray spectrum directly reveals that the soft excess 
is not composed of narrow emission lines and that it has significant spectral curvature. A power-law model fails 
to fit it, while a blackbody produces a nearly acceptable fit. We compare the observed spectral variability with the 
results of previous studies of NGC 405 1 . 

Subject headings: galaxies: active — galaxies: nuclei — galaxies: individual (NGC 4051) — galaxies: Seyfert 
— X-rays: galaxies — ultraviolet: galaxies 



1. INTRODUCTION 

NGC 4051 is a Narrow-Line Seyfert 1 (NLS1) galaxy with 
a rich observational history. Its optical emission-line spec- 
trum was first noted by Hubble (1926), and it is one of the 
archetypical Seyfert galaxies (Seyfert 1943). It has been stud- 
ied with every X-ray mission since Einstein, as well as with 
EUVE, HST and WE, because it is bright [V w 13.5; z = 
0.002295 ±0.000043 based on optical emission lines (de Vau- 
couleurs et al. 1991)] and has relatively little Galactic absorp- 
tion [N H = (1.3 ±0.1) x 10 20 cm" 2 ; Elvis, Lockman, & Wilkes 
1989]. NGC 4051 is highly variable in X-rays and has even 
been seen virtually to shut off (decreasing in flux by a factor of 
k 20 from the average; e.g., Guainazzi et al. 1998; Uttley et al. 
1999). Its typical 2-10 keV luminosity is (2-5)xl0 41 erg s" 1 
(Hq = 70 km s" 1 Mpc" 1 ; qo = 0.5). At lower energies the X-ray 
continuum is dominated by a variable soft excess (e.g., Turner 
& Pounds 1989; Guainazzi et al. 1996, hereafter G96). 

NGC 405 1 is also known to contain a 'warm absorber' (X-ray 
absorbing highly ionized gas) from ROSAT and ASCA observa- 
tions (e.g., McHardy et al. 1995; G96; Komossa & Fink 1997; 
Reynolds 1997; George et al. 1998). Although edges from 

VII and O VIII have been statistically detected in many re- 
cent analyses of NGC 4051, the addition of these features can- 
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not entirely account for the observed spectral complexity (even 
when a blackbody model for the soft excess is also included 
in the fitting). For example, G96 identified a narrow spectral 
feature near 1 keV but were unable to determine if this feature 
was an absorption edge or emission line. George et al. (1998) 
suggested that the strong emission lines that are expected to 
arise in the warm absorber (e.g., Krolik & Kriss 1995; Netzer 

1996) may comprise a significant fraction of the observed soft 
excess and contribute to the spectral complexity; however, all 
X-ray observations to date have lacked the spectral resolution 
needed to confirm or reject this hypothesis. G96 also argued 
for significant variability of the O VII edge while the O VIII 
edge remained roughly constant in strength, a somewhat sur- 
prising result if the O VIII absorption occurs closer to the cen- 
tral source on average than the O VII absorption (see Reynolds 

1997) . The best-fit ionization parameter of the warm absorber 
in NGC 4051 does not always track source intensity; such be- 
havior could arise if the warm absorber is not in photoioniza- 
tion equilibrium or has a multi-zone nature (e.g., McHardy et al. 
1995; Nicastro et al. 1999). Finally, we note that the possible 
presence of intrinsic UV absorption was briefly suggested by 
Voit, Shull, & Begelman (1987), but to date such absorption 
has not been studied in any detail. 
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In this paper, we present the results from our analyses of si- 
multaneous Chandra High Energy Transmission Grating Spec- 
trometer (HETGS) and HST Space Telescope Imaging Spectro- 
graph (STIS) observations of NGC4051. Through this study, 
we aim to constrain the dynamics (e.g., bulk velocity and veloc- 
ity dispersion) and geometry (e.g., location and covering factor) 
of the X-ray and UV absorbers, and to determine whether these 
features might plausibly arise in the same gas (e.g., Mathur, 
Elvis, & Wilkes 1995). As our observation comprises the first 
X-ray grating spectroscopy of NGC 405 1 , we shall also use the 
high spectral resolution of the HETGS to place the best con- 
straints to date on the spectral form and nature of the soft excess 
and to examine flux and spectral variability. In § 2 we describe 
the observations and data analysis, and in § 3 we interpret our 
results. 

2. OBSERVATIONS AND DATA ANALYSIS 

2. 1 . Chandra HETGS Observation 

2.1.1. Observation Details and Basic Analysis 

NGC 405 1 was observed by Chandra (Weisskopf et al. 2000) 
on 2000 April 24-25 (starting at 06:08:25 UT) using the 
HETGS (C.R. Canizares et al. 2001, in preparation). 9 This 
observation was part of the Chandra Cycle 1 guest observer 
program (observation identification number 00859). The detec- 
tor was the Advanced CCD Imaging Spectrometer (ACIS; GP. 
Garmire et al. 2001, in preparation). 9 The observation was con- 
tinuous with a total integration time of 81.5 ks. We used the 
Chandra Interactive Analysis of Observations (CIAO) software 
Version 1.1.4 (M. Elvis et al. 2001, in preparation) to reduce 
the data. 10 

The zeroth order X-ray spectrum shows substantial photon 
pileup. The nuclear position derived from the Chandra image is 
a 2 ooo = 12 h 03 m 09?58, <5 20 oo = +44 d 31 m 52?9, in agreement with 
the radio position from Ulvestad & Wilson (1984); our position 
is offset from theirs by 0."6 to the North, while the absolute 
Chandra astrometry for our observation is expected to be good 
to within 0."6. n Comparing the point spread function (PSF) 
wings of the NGC 4051 image to the PSF of a point source 
(the HETGS observation of Capella), we find good agreement 
with no clear evidence for an extended circumnuclear X-ray 
component in NGC 4051. This is in contrast to the findings 
of Singh (1999) who claimed to detect extended circumnuclear 
emission in a ROSAT HRI observation; the extended emission 
was claimed to contain 21 ± 6% of the total flux. Our obser- 
vation has > 5 times better spatial resolution than the ROSAT 
HRI observation and is « 8 times longer (the ROSAT HRI in- 
strument had a comparable effective area to the zeroth order of 
the HETGS). The extended emission reported by Singh (1999) 
may be a result of limitations in the ROSAT 'de-wobbling' pro- 
cedure or the theoretical model of the ROSAT HRI PSF. We 
also comment on the presence of two weak off-nuclear (« 17" 
and « 43") sources, located at a 20 oo = 12 h 03 m ll?04, <5 20 no = 
+44 d 31 m 46. s and a 200 o = 12 h 03 m 13?51, S 20 qo = +44 d 31 m 46?7, 
coincident with the galactic disk of NGC 405 1 . These sources, 
which do not have any apparent discrete counterparts in Palo- 
mar Observatory Sky Survey (POSS) or ffi>r images, have w 24 
and w 1 1 counts, respectively, corresponding to rest-frame lu- 

9 For additional information on the HETGS and ACIS, see the Chandra Proposers' 

10 See http://asc.harvard.edu/ciao. 

11 See http://asc.harvard.edu/mta/ASPECT. 

12 See http://space.mit.edu/ASC/calib/hetgcal.html. 



minosities of Lo.5_g kev ~ 1.4 x 10 38 erg s 1 and Lo.5-8 kev ~ 
6.6 x 10 37 erg s" 1 if they are at the distance of NGC 4051. They 
are plausibly X-ray binaries or luminous supernova remnants in 
NGC 4051, although they could also be unrelated background 
sources; the probability of obtaining one background source co- 
incident with the optical disk of NGC 405 1 is w 40% at the ob- 
served flux levels. 

The HETGS produces higher order spectra from two grat- 
ing assemblies, the High Energy Grating (HEG) and Medium 
Energy Grating (MEG); each produces positive and negative 
orders. The two first-order MEG spectra agree well with each 
other and with the first-order HEG spectra. The averaged first- 
order MEG and HEG spectra have signal-to-noise ratios (S/Ns) 
of s=s 4.6 and w 2.3 (near 7 A) for bin sizes of 0.005 A and 
0.0025 A, respectively. The higher orders have only a small 
number of counts, and we therefore exclude them from our 
analysis. We have used CIAO to produce Ancillary Response 
Files (ARFs) in order to flux calibrate the spectra, and we have 
corrected for the Galactic absorption and the cosmological red- 
shift. The flux calibration is currently estimated to be accu- 
rate to w 30%, 20% and 10% in the 0.5-0.8, 0.8-1.5 and 1.5- 
6 keV bands, respectively (H. L. Marshall 2000, private com- 
munication). 12 We have also corrected for the wavelength shift 
of 0.054% induced by the thermal contraction of ACIS-S pix- 
els (H. L. Marshall 2000, private communication). We present 
the mean MEG first-order spectrum in Figure 1 . This spectrum 
contains a total of 32087 counts between 1 .8 A and 25 A (0.5- 
6.9 keV). The MEG velocity resolution is 1400 kms" 1 at 5 A 
and 350 kms" 1 at 20 A. 

The general shape and flux of the X-ray spectrum are con- 
sistent with previous observations of NGC 4051. The aver- 
age fluxes measured from the Chandra spectrum are (1.4 ± 
0.3)xl0" n erg cm" 2 s" 1 in the 0.5-2 keV band and (1.2 ± 
0.1)xlO -11 erg cm" 2 s" 1 in the 2-7 keV band; the errors on 
these figures represent the flux calibration uncertainties men- 
tioned above. We also observe short-timescale variability (see 
Figure 2) and a state of low flux during the last 15 ks of the 
observation. For this interval, the average Chandra fluxes are 
(2.9±0.6)xl0 -12 erg cm" 2 s" 1 in the 0.5-2 keV band and 
(5.5 ±0.6)xl0~ 12 erg cm" 2 s" 1 in the 2-7 keV band; these 
fluxes, while quite low for NGC 4051, are still about four times 
higher than during the most extreme of the ultra-dim X-ray 
states mentioned in § 1. The largest amplitude of count-rate 
variability observed is a factor of 12.6 ±0.7, and the average 
count rate during the last 15 ks of the observation is 4.8 ±0.1 
times smaller than for the preceding data. We note that the 
low-flux state could not have lasted for more than « 3 days, 
as NGC 405 1 was observed to be in a more typical flux state 
shortly after our observation (P. Uttley 2000, private communi- 
cation). 

2.1.2. Detected X-ray Spectral Features 

The most interesting features revealed in the high-resolution 
X-ray spectrum of NGC 405 1 are narrow absorption and emis- 
sion lines, similar to those recently detected in other Seyfert Is 
such as NGC 3783 (Kaspi et al. 2000) and NGC 5548 (Kaastra 
et al. 2000). Because only a small number of the individual 

Observatory Guide at http://asc.harvard.edu/udocs/docs. 
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FIG. 1. — Chandra MEG first-order spectrum of NGC 4051. The bin size is 0.005 A, and the spectrum has been smoothed using a boxcar filter that is three bins 
in width. Fluxes are corrected for Galactic absorption. Dots with error bars show the typical ±1ct statistical error at various wavelengths. The identified absorption 
(a) and emission (e) features are marked. Strong unidentified features are indicated. Some additional spectral complexity appears to be present, particularly in the 
middle panel (e.g., 14-16 A); however, the wavelengths of the apparent features do not seem to correspond to strong expected absorption or emission features from 
either of the two identified X-ray absorption systems. This additional complexity may be the cumulative result of a forest of individually weak lines (e.g., Kaspi 
et al. 2001). Note the difference in vertical scale between the top two panels and the bottom panel. 




Time [ks] 

FIG. 2. — Chandra MEG light curve and hardness ratio for NGC 4051. The 
hardness ratio is defined to be the ratio of the 2-8 keV and 0.5-2 keV counts. 
The top panel has a bin size of 300 s, and the bottom panel has a bin size of 
3000 s. Note the low count rate during the last rss 15 ks; the spectrum clearly 
hardens during this time. The horizontal lines in the top panel show the times 
when STIS data were acquired. 



absorption features are statistically significant at a high level of 
confidence, we have adopted the following procedure to mea- 
sure relevant properties of the lines. We created 'velocity spec- 
tra' by adding, in velocity space, several absorption lines from 
the same ion. This method has the advantage of a significant 
increase in S/N and the disadvantage of somewhat deteriorated 
spectral resolution since, for each ion, the resolution is then 
defined by the line with the shortest wavelength. The veloc- 
ity spectra were built up on a photon-by-photon basis from the 
raw data, rather than by interpolating spectra already binned in 
wavelength. In all cases, the two strongest predicted features 
for each ion were included (H. Netzer 2000, private communi- 
cation), and in some cases other lines were added in order to 
improve the strength of the signal. We chose the lines to be rel- 
atively free from contamination by unrelated adjacent features 
based upon observations and modeling of other Seyfert galax- 
ies with warm absorbers (mainly NGC 3783; Kaspi et al. 2001). 
In addition to combining the absorption lines from single ions, 
we also combined all the transitions from H-like ions into one 
spectrum and all the transitions from He-like ions into another. 
This step was phenomenologically motivated by the observed 
similarities between many of the spectra for the H-like and He- 
like ions. The velocity spectra created in this way are shown 
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FIG. 3. — Chandra MEG velocity spectra showing coadded lines from H-like and He-like ions of O, Ne, Mg and Si. The bin size is 400 km s -1 . The wavelengths 
of the lines coadded in each spectrum are listed in the respective panels, except for the total H-like and He-like spectra which contain all of the lines listed for that 
type of ion. Absorption systems at two distinct velocities (-2340 km s -1 and -600 km s _1 ) are detected when all ions are considered (compare the upper two panels). 



in Figure 3. We measured the velocity shifts and widths of 
the observed absorption systems from Gaussian fits to the total 
H-like and He-like spectra. We detect two distinct absorption 
components with > 99.9% confidence according to the F-test 
(see Bevington & Robinson 1992): a high-velocity system at 
-2340 ± 130 km s" 1 (Ax 2 = 59.70 in the H-like spectrum) and 
a low-velocity system at -600 ± 130km s" 1 (A^ 2 = 38.34 in the 
He-like spectrum and Ax 2 = 15.41 in the H-like spectrum). The 
two absorption components have FWHMs of 980 ± 3 10 km s" 1 
and 770 ±460 km s" 1 , respectively. These values are consistent 
with the resolution of the MEG at the wavelengths of some of 
the lines in each coadded spectrum. 

Constraining the velocity shifts of the two components to 
have the values we measured from the H-like (high-velocity 
component) and He-like (low-velocity component) ions, we 
then performed F-tests on the velocity spectra of the individ- 
ual ions in order to check statistically for the presence of both 



absorption systems from each ion. This process consisted of 
fitting the continuum from -5000 to 5000 km s -1 with a linear 
model and noting the change in x 2 when narrow (unresolved) 
Gaussian absorption lines were added at fixed velocities cor- 
responding to either high-velocity or low-velocity absorption. 
Because the only free parameter in the Gaussian model was 
the normalization, the values of Ax 2 correspond to one pa- 
rameter of interest. The results of this procedure are shown 
in Table 1. According to the F-test, we detect absorption at 
the > 99.9% confidence level from the high- velocity system in 
Si XIV, Mg XII, and Ne X, and from the low-velocity system in 
O VII. Absorption is also suggested from the remaining H-like 
or He-like species of O, Ne, Mg and Si. Because deviations in 
the continuum from a linear model may introduce small inac- 
curacies into our F-test results, we only consider features indi- 
cated at the > 99.9% confidence level to be reliable detections. 
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TABLE 1 

Values of Ax 2 for Co-added X-ray Line Profiles 





Low Vel. 


High Vei. 


Ion 


A X 2 


Ax 2 


Si XIII 


4.66 


0.38 


Si XIV 


0.49 


14.75 


MgXI 


4.72 




MgXII 


0.38 


18.05 


NelX 


5.82 


5.77 


NeX 


3.44 


21.19 


O VII 


37.72 


0.49* 


ovm 


6.79 


3.44* 



NOTE. — Values of Ax 2 represent the difference in x 2 derived from 
linear models and linear models with one Gaussian component to rep- 
resent absorption from the indicated velocity system, except for the val- 
ues marked with '*'. In these cases, Ax 2 is the difference in x 2 de- 
rived from linear models with one Gaussian component to represent low- 
velocity absorption and linear models with Gaussian components to rep- 
resent absorption at both velocities. All fits have at least 40 degrees of 
freedom. For 40 degrees of freedom, some Ax 2 values and their corre- 
sponding significance levels are: Ax 2 = 4.08 / 95%; Ax 2 = 7.31 / 99%; 
Ax 2 = 12.6 / 99.9% (see Table C.5 of Bevington & Robinson 1992). 

In Table 2 we list information about the individual absorption 
and emission lines. We have used the same procedure detailed 
in the preceding paragraph to search for individual absorption 
lines in the X-ray spectrum, again constraining the absorption 
components to be narrow and to occur at wavelengths corre- 
sponding to the velocities given above. Although this procedure 
was very effective at identifying strong features, we note that 
it cannot entirely account for the observed spectral complexity 
(see Figure 1 for clarification). In regions of the spectrum where 
there are not enough counts to permit x 2 fitting, we have man- 
ually searched for absorption features. Similarly, we have not 
performed x 2 statistical significance testing on any of the iden- 
tified emission lines because many of these occur in regions of 
the spectrum with low S/N. We measured the velocity shifts and 
widths of the emission lines from simple Gaussian fits. Where 
allowed by the strength of the signal, we have calculated the er- 
rors on these values at the 90% confidence level (A% 2 = 2.71). 
In other cases, we have estimated the errors by visual inspec- 
tion of the features. We note that the emission lines appear to 
be unresolved (see § 3.1 for further discussion) and are gen- 
erally consistent with zero velocity within the NGC4051 rest 
frame. We present the X-ray emission-line fluxes in Table 3. 
In general, the statistical significances of the emission lines and 
the low S/N absorption lines can be judged by the sizes of the 
errors on the equivalent widths (EWs). In order to measure the 
EWs of all the individual features, we have excised the line- 
dominated areas in the spectrum and fit the remainder with a 
cubic spline which serves as the continuum reference point in 
the centers of the lines. Errors on EWs represent the combina- 
tion of a global 10% uncertainty in the continuum level and the 
statistical photon noise. 

Some of the strong absorption lines occurring at the long- 
wavelength end of the MEG spectrum appear to be resolved 
(see Figure 4). Fitting a Gaussian profile to the coadded 
O VII velocity spectrum yields a best-fit FWHM for the low- 
velocity absorption system of 850^0 km s" 1 (errors are for 
Ax 2 = 2.71), whereas the instrumental resolution (FWHM) of 
the MEG at 17.768 A (the line in the O VII velocity spectrum 
with the shortest wavelength and therefore the worst resolution) 
is 390 km s -1 . Correcting for the instrumental broadening gives 
a true FWHM of 750!$$ km s" 1 for the low-velocity X-ray ab- 
sorption system. Although the high-velocity system appears to 



TABLE 2 

X-ray Absorption and Emission Lines From NGC 405 1 



Ion 


Rest 
A 

(A) 


Ax 2 


Velocity 
(System 
or km s -1 ) 


Measured 
EW 
(mA) 


Absorption Lines 


o A. V 1 


A 798 
4. / Zo 


z.zu 


Hign 


7 r\+^- z 

- 5 A 


Si Al V 


6.181 


13.05 


High 


13.0";* 


Si XIII 


6.648 


3.58 


Low 


5-3111 


Mg XII 11 


8.419 


14.14 


High 


18-71J! 


NeX a 


10.239 


2.90 


High 


6 2+ 2 - 4 


NeX a 


12.132 


4.60 


Low 


1 9+ 4 - 1 


NeX a 


12.132 


11.52 


High 


35 9 +6 ' 5 


Fe XVIP 


12.263 


7.30 


Low 


16.2 ±4.0 


NelX 


13.447 


4.70 


Low 


15.9+5 i 


NelX 


13.447 


2.90 


High 


n.o+i? 


Fe XVIII 


14.539 


2.33 


Low 


9-9% 


O VHP 


16.006 


3.48 


Low 


10.3 ±3.2 


O VHP 


16.006 


2.93 


High 


6.3 ±3.6 


O VII 


17.768 


3.15 


Low 


17.5 ±8.5 


O VII 


18.628 


6.28 


Low 


24.8 ±9.5 


O VIII 


18.967 




Low 


30.1 ±4.9 


O VII 


21.602 




Low 


37.1 ±12.2 


Emission Lines 


Fe I -XVIII 


1.937 




-250 ±750" 


47 7+13.4 


Si XIII 


6.741 




q+250 


12.0_f4 


NeX 
NelX 


12.132 
13.553 




_ 5 >o° 


15 l +6: ' 
9-7-f? 


NelX 
O VIII 
O VII 
O VII 


13.699 
18.967 
21.807 
22.102 




180f| 
_ 6() +ioo 

-122+'^ 


48.6_gf 
21 4+" 
37.3+lJI 
180.l2|| 



"Possibly blended with Fe absorption lines. 

b The velocity of Fe K« is not well-constrained due to the low S/N 
in that region of the spectrum (see Figure 12). The velocity reported 
here is the average velocity measured in the MEG and HEG spectra. 
The wavelength of this observed line corresponds to the expected 
iron Ka line from Fe I to Fe XVIII. 



TABLE 3 

X-ray and UV Emission-Line Fluxes 





Rest 




Ion 


Wavelength 

(A) (xur 


Flux 

44 erg cirr 2 s~') 


X-ray Emission Lines 


Fe I -XVIII 


1.937 


Jz -10 


Si XIII 

NeX 


6.741 
12.132 


1 3+ 07 
1.0 ±0.4 


NelX 


13.553 


0.6 ±0.4 


NelX 
OVIII 


13.699 
18.967 


? q+0.9 
z -*-0.8 

1.6 ±0.7 


O VII 


21.807 


2.4 ±1.4 


O VII 


22.102 


11 ± 3 


UV Emission Lines 


HI 


1215.7 


137 a 


NV 


1238.8/1242.8 


24" 


Si II 


1304.4 


12±3 


CII 


1334.5 


3.1 ±0.5 


Si IV 


1393.8/1402.8 


21 ±5 


NIII 


1485.8 


3.0 ±0.7 


CIV 


1548.2/1550.8 


117" 


He II 


1639.8 


18±5 



"Fluxes for the heavily absorbed UV lines, Lya, N V, and C IV, 
are highly uncertain as a result of the strong absorption. 



be unresolved in all of the coadded line profiles, we can use the 
Ne X profile to place an upper limit on its measured FWHM of 
1 140-1860 km s" 1 , which corresponds to a true width of <910- 
1730 km s" 1 . We are unable to use the HEG spectrum to place 
tighter constraints on the lines' widths due to its low S/N. 
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FIG. 4. — O VII velocity spectrum from Figure 3 built up from the 17.768 A, 
18.628 A and 21.602 A lines. The top panel shows a simulated unresolved ab- 
sorption line profile as it would appear in the MEG at a wavelength of 17.768 A 
(the shortest wavelength O VII line and therefore the one with the worst MEG 
velocity resolution). The coadded spectrum is shown with two different bin 
sizes to allow easy comparison with the instrumental response. 

We have checked the two strongest emission lines, O VII 
22.102 A and Ne IX 13.699 A, for flux variability between the 
high and low states (see § 2.1.3) and find no clear evidence for 
it. The upper limit on the variability amplitude of the O VII line, 
which is stronger than the Ne IX line, is w 51%. 

2.1.3. Spectral Variability and Soft X-ray Excess 

We have independently analyzed the grating spectra from the 
first w 65 ks of the observation (hereafter the 'high state') and 
the last « 15 ks ('low state') in order to study spectral vari- 
ability (see Figure 2). We have binned the spectra to have a 
minimum of 15 counts per data point, in order to allow x 2 
fitting. Since we have flux calibrated the spectra and cor- 
rected for Galactic absorption and redshift, we have used the 
CIAO tool SHERPA to fit spectral models directly to the data, 
without having to make further accounting for instrumental re- 
sponse. We simultaneously fit all four first-order HETGS spec- 
tra, constraining all shared spectral shape parameters to be the 



same, but allowing the absolute model normalizations to vary 
(to allow for absolute instrumental flux-calibration uncertain- 
ties). We note that the model normalizations were in general 
consistent to within 10%. In order to obtain the most reliable 
constraints on continuum spectral properties, we have excised 
all strong narrow features prior to fitting. In the high-state and 
average spectra, we have removed all the lines identified in Ta- 
ble 2 as well as the unidentified features marked in Figure 1 . In 
the low state, we have removed only those parts of the spectrum 
containing strong emission lines. 

We fit the > 2 keV continuum with a power law. The fits 
were statistically acceptable (see Table 4), and the best-fit pho- 
ton indices with 90% confidence ranges (Ax 2 = 2.71) were 
1.85 ±0.05 for the high state and 0.83ig;}| for the low state. 
The flattening of the hard power-law continuum in the low state 
is qualitatively consistent with the behavior seen in previous 
X-ray observations (see §3.2 for discussion). Upon extrapo- 
lating the power-law fits to lower energies, the soft excess be- 
comes strongly apparent (see Figure 5). It appears to be con- 
tinuous rather than composed of many narrow emission lines. 
Given the detection of this strong and continuous soft excess in 
high-resolution X-ray data, it seems unlikely that the extreme- 
ultraviolet to hard X-ray continuum of NGC 405 1 is composed 
only of a single power law (see Uttley et al. 2000). Compari- 
son of the two panels of Figure 5 shows that the soft excess is 
rapidly variable. In order to measure the amplitude of this vari- 
ability, we have subtracted the power-law models shown in Fig- 
ure 5 from the data and calculated the remaining 0.5-1.5 keV 
fluxes in both the high and low states. The ratio of the high- 
state to low-state soft excess fluxes we obtain is 2.56+q g!?. The 
soft excess is proportionally stronger (relative to the underlying 
continuum) in the low state than in the high state; the fraction 
of the 0.5-1.5 keV flux due to the soft excess is « 76% in the 
low state and w 36% in the high state. We comment that our ba- 
sic results on the soft excess are not sensitive to the correction 
for Galactic absorption; NGC 4051 has a small and precisely 
measured Galactic column density (see §1). 

We first attempted to fit the soft excess with a power law. 
This model is rejected at > 99.9% confidence for the high state 
and at 98.7% confidence for the low state (see Table 4). In the 
high state, this fit displays systematic residuals in the soft band, 
indicating that the soft excess has significant spectral curvature. 
Therefore we next attempted to fit it with a curved model, a 
blackbody. This provides a highly significant improvement to 
the fit in the high state (Ax 2 = 117.7), while in the low state 
the fit improves only slightly (Ax 2 = 3.8); the shape of the soft 
excess is poorly constrained in the low state due to its low S/N. 
The addition of a second blackbody component with a differ- 
ent temperature does not result in a statistically significant im- 
provement in either state. Although the single blackbody model 
is still statistically rejected with high confidence, this is largely 
the result of unmodeled spectral complexity, much of which 
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X 2 /d.o.f. (Degrees of Freedom) for X-ray Spectral Fitting 





Power law 


Power law and 


Power law and 




(> 2 keV data only) 


power law 


blackbody 


Low State 
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High State 
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Average Spectrum 
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2133.7/2019 
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FIG. 5. — (a) High-state and (b) low-state MEG negative first-order spectra of NGC4051. Each data point represents a minimum of 60 counts in the top panel 
and a minimum of 15 counts in the bottom panel. Strong absorption and emission lines have been excised from the spectra. We have fit the continua above 2 keV 
with power laws (using all four first-order HETGS spectra; see §2.1.3), which we have extrapolated to lower energies to show the soft excess. Note the significant 
spectral curvature of the soft excess apparent in the top panel (at < 0.7 keV). 



may be from iron L-shell lines that are individually below the 
detection threshold but which have a significant cumulative ef- 
fect (e.g., Kaspi et al. 2001). The spectral shape of the black- 
body component is consistent to first order with the shape of 
the soft excess. The best-fit blackbody temperatures (with 
90% confidence uncertainties) for the high and low states are 
0.1045 ±0.0006 keV and 0.1 126 ±0.0019 keV, respectively. 
Adding an intrinsic neutral absorption column does not improve 
the fit; the upper limit on the column density of such a compo- 
nent is w 10 20 cm" 2 . We note that the presence of any intrin- 
sic, neutral absorption in NGC4051 would only increase the 
strength of the required soft X-ray excess. 

Adopting the blackbody as the best model for the soft-excess 
continuum, we then statistically tested for the presence of ab- 
sorption edges in the high-state spectrum. We find evidence 
for an O VII absorption edge (A% 2 = 25.5) with a velocity con- 
sistent with the range of velocities encompassed by the X-ray 
absorption lines and with r = 0.220 ± 0.035, implying an O VII 
column density of (8 ± 1) x 10 17 cm -2 . The precise energy of 
the edge (0.741 ± 0.019 keV) is not well constrained due to the 
loss of resolution necessary to reach x 2 statistics. Our test re- 
sults also suggest the presence of an O VIII edge (A% 2 = 17.1) 



with a similar and equally uncertain implied velocity (energy 
in the range 0.8714-0.8787 keV) and with r = 0.157 ±0.032, 
implying an O VIII column density of (1.4 ± 0.3) x 10 18 cm -2 . 
In addition, there is statistical evidence for another absorption 
edge with an energy of « 0.91 keV and r « 0.17 (A% 2 ss 20; the 
energy, strength and significance depend on the presence or ab- 
sence of other edges in the model). This feature may be caused 
by the previously mentioned iron L-shell absorption lines. We 
note that the presence of such absorption lines may also result 
in systematic uncertainties in the O VII and O VIII edge depths 
and hence the implied column densities. 

2.2. HST STIS Observation 

2.2. 1 . Observation Details and Basic Analysis 

NGC4051 was observed with STIS (Woodgate et al. 1998) 
on HST during the Chandra observation (see Figure 2). This 
observation was part of the HST Cycle 8 guest observer pro- 
gram (data sets O5F001010, O5F001020, O5F001030, and 
O5F001040). The STIS observation was made with the 
medium-resolution echelle (E140M mode), using an aperture 
size of 0."2 x 0."2. The total exposure time was 10.3 ks (4 or- 
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FIG. 6. — Heavily binned STIS spectrum of NGC4051. The bin size is 0.25 A. Intrinsic features are labeled above the spectrum. Strong Galactic features are 
labeled below the spectrum (wavelengths are from Morton, York, & Jenkins 1988). 



bits), and the wavelength coverage was 1173-1730 A with a 
resolution of « 7 km s" 1 and an absolute wavelength calibra- 
tion good to within one resolution element. We have analyzed 
the standard STScI pipeline spectrum and compared it with the 
spectrum produced using the Interactive Data Language (IDL) 
software of the Instrument Definition Team (IDT) for STIS 
(see, e.g., Crenshaw et al. 2000). We find good overall con- 
sistency between the two spectra, and we use the IDT spec- 
trum in the analysis below. Figure 6 shows the UV spectrum of 
NGC4051 from 1190-1720 A. 

We have searched for orbit-to-orbit variability of the UV con- 
tinuum, absorption lines and emission lines. We detect none 
despite the simultaneous large-amplitude X-ray variability (see 
Figure 2; compare with Done et al. 1990). The continuum 
fluxes in broad wavelength bands are consistent to within « 3%. 
This lack of detected variability justifies our use of the mean 
UV spectrum below. 

2.2.2. Detected Ultraviolet Absorption Lines 

The STIS spectrum of NGC 405 1 reveals the presence of a 
number of intrinsic absorption systems with velocities ranging 
from -650 km s" 1 to 30 km s" 1 . The velocity widths (FWHM) 
of the components range from 23 km s" 1 to 165 km s" 1 . We 
see up to nine major systems in C IV and N V, and we detect 
some of these in Si IV, Si III, Si II, and C II as well (see Table 5 
and Figures 6, 7, and 8). We note, however, that our criteria for 
identifying distinct absorption systems are somewhat subjec- 
tive. We classified two adjacent absorption systems as distinct 
if they appeared visually separate in at least one of the tran- 
sitions shown in Figure 7. Unfortunately, we were unable to 
examine intrinsic Lya absorption from NGC 4051, since it dis- 
plays a saturated absorption trough due to Galactic and intrinsic 
absorption (as well as perhaps intergalactic material) that ex- 
tends from about -980 km s" 1 to 50 km s" 1 (in the rest frame of 
NGC 405 1 ; see Figure 9). This trough includes the entire range 



of observed intrinsic UV absorption velocities and is therefore 
consistent with our other results. We note that there are also 
two Lya absorption systems that are not saturated with veloc- 
ities of ss 110 km s" 1 and w 260 km s" 1 . These systems do 
not seem to have counterparts in any of the high-ionization UV 
lines; they may represent low-ionization, high-velocity clouds 
in NGC 4051 (as in our Galaxy; e.g., Wakker & van Woer- 
den 1997). We identify a Galactic absorption system candidate 
(having velocity approximately consistent with Galactic origin) 
in all ions but N V (labeled 'G' in Figures 7 and 8), and we note 
that this absorption could possibly be subdivided into two dis- 
tinct systems. We also point out that system 1 could be Galactic 
in origin. 

In Figure 7 and Table 5 we present the absorption systems 
from the higher ionization ions (C IV, N V, and Si IV). In order 
to measure the velocity shifts, FWHMs, and EWs of the sys- 
tems (many of which are saturated), we have fit the continuum 
on either side of the absorption with a cubic spline. The rele- 
vant continuum fits are shown in Figure 7. Although this fitting 
technique may underestimate the continuum level if there is a 
significant narrow component to the emission, visual inspection 
of the less-strongly absorbed emission lines suggests that this 
is probably not a highly significant effect. Because the error on 
the spline-fit continuum level is difficult to quantify, we adopt 
a global uncertainty of 10% on these values, which is consis- 
tent with the range of possible continuum fits we considered. 
The errors on the equivalent widths represent a combination of 
this uncertainty in the continuum level and the actual photon 
noise. In general, where two or more absorption components 
appear to be strongly blended together, we report the sum of 
the EWs rather than attempting to separate them. The veloci- 
ties and velocity FWHMs of the systems have been determined 
from simple Gaussian fits to the absorption. The velocity shifts 
reported in Table 5 represent the averages from N V, C IV, and 
Si IV; we take the root-mean-square deviations of these mea- 
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FIG. 7. — Intrinsic UV absorption systems in NGC4051. The bin sizes for the upper, middle, and lower panels, respectively, are 0.016 A, 0.013 A, and 0.015 A. 
Spectra have been smoothed with a boxcar filter that is five bins in width. The absorption features are marked 'G' for Galactic or numbered 1-9 for the intrinsic 
systems. The cubic spline continuum fits are shown. Note that many of the absorption lines are saturated. 



surements (which are typically 10-15 km s" 1 ) to be the best 
estimates of the uncertainties involved. 

In Figure 8 we present the absorption systems from the 
lower-ionization ions (C II, Si III, and Si II). In these lines we 
identify an additional absorption system (labeled as '10') which 
has a velocity of « -80 km s" 1 . For the lower-ionization ions, 
we tentatively identify counterparts to some of the absorption 
systems found using the higher-ionization ions (the statistical 
significances of these features are strongly dependent upon the 
adopted continuum level). The most notable counterparts are 
systems 1, 4, 7, 8, and 9, which are found in all ions, sys- 
tem 2 (found in Si II and Si III), and system 5 (found in C II 



and Si III). The low-velocity absorption systems (8-10) may 
originate within the host galaxy rather than the active nucleus. 
However, the high-velocity absorption systems (2-5) seen in 
the low-ionization ions are most likely related to the nucleus. 
Such high-velocity absorption systems in low-ionization ions 
have been observed in only one other Seyfert galaxy, NGC 415 1 
(Weymann et al. 1997). 

2.2.3. Comparison with IUE and Ultraviolet Emission Lines 

Peterson et al. (2000) present He II and C IV emission lines 
for NGC 405 1 from a mean UV spectrum of 3 1 observations 
carried out by the International Ultraviolet Explorer (WE) 
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FIG. 8. — Intrinsic UV absorption systems in NGC4051 from the low- 
ionization ions C II, Si III, and Si II. The bin size is f» 0.15 A. Spectra have 
been smoothed with a boxcar filter that is five bins in width. The expected 
wavelengths of the absorption systems from Figure 7 and Table 5 are marked 
as vertical lines in each panel. Some detected absorption features are labeled 
as in Figure 7. Here we also identify an additional absorption system, labeled 
as '10', which has a blueshifted velocity of -80 km s . 



15 - 



■g 10 
o 

60 
1) 



2 5 
X 

I 




1209 



i i i i i 

1215 

Rest Wavelength (A) 



t t 

I i I i I i I i I i I i I i I i I i I i I 
1212 1215 1218 1221 



1224 



FIG. 9. — Absorption near Lya in the spectrum of NGC4051. The bin size 
is 0.13 A. The spectrum has been smoothed with a boxcar filter that is five bins 
in width. Note the wide absorption trough and the narrow geocoronal emission 
line near 1213 A (the spectrum has been corrected for the cosmological red- 
shift). The two absorption systems with redshifted velocities of fs 1 10 km s _1 
and Kb 260 km s _1 are marked with arrows. These systems may originate in 
low-ionization, high- velocity clouds within NGC 405 1 (as in our Galaxy). 




I I I I I i 

1500 1520 1540 1560 1580 1600 

Observed Wavelength [A] 

FIG. 10. — Mean IUE spectra of NGC 4051 (around the C IV emission line) 
for several epochs. The epochs are presented on the left, and the number of 
spectra averaged is on the right. Note the absorption fa 2.5 A blueward of 
the line peak, which is present in all epochs and is consistent with the absorp- 
tion seen in the HST STIS spectrum. The small wavelength shifts between the 
spectra are probably due to wavelength calibration problems and the low res- 
olution. These small wavelength shifts are one reason that the absorption is 
difficult to see in the mean spectrum presented by Peterson et al. (2000). 

Short-Wavelength Prime camera (SWP; wavelength range 
~ 1150-2000 A) over a period of 16 years. Though it is not 
clear from Peterson et al. (2000), a careful plotting of the C IV 
line reveals the likely existence of absorption in 90% of the 
spectra (see Figures 10 and 11; the lack of apparent absorption 
in the other spectra is probably due to their low S/N). Thus it 
appears that C IV absorption has existed in NGC 405 1 since at 
least 1978. 

The STIS spectrum of NGC 405 1 is consistent with the mean 
IUE spectrum at wavelengths longer than 1225 A. Around the 
Lya line the IUE spectrum is contaminated by geocoronal Lya 
emission, and no comparison is possible. The emission lines 
from both spectra agree well in their profiles. The blueward 
asymmetry of the He II and C IV broad emission described by 
Peterson et al. (2000) is also detectable in the STIS spectrum. 
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FIG. 1 1 . — Absorption and emission systems as seen by (a) and (b) STIS, (c) IUE (mean of the 1978-1988 and 1994 epochs presented in Figure 10), and (d) 
Chandra; panel (b) represents the STIS C IV line binned to the velocity resolution of the MEG. The approximate velocity resolution is given for each plot. The 
velocity units in the panels showing C IV correspond to the 1550.8 A line, and dots with error bars mark the velocities of the two X-ray absorption components 
resolved by Chandra. Comparison of the panels stresses the point that the X-ray warm absorber may be subdivided into further systems that cannot be resolved with 
current X-ray instruments. Note also that no strong ultraviolet absorption is seen corresponding to the high-velocity X-ray absorption component. 



This effect is also seen in other NLSls (e.g., Leighly 2000). 
We list the UV emission-line fluxes in Table 3; in all cases the 
fluxes have been measured from unabsorbed line profiles fit to 
the data as described in § 2.2.2. Note that the fluxes given for 
Lya, N V, and C IV are quite uncertain due to their heavily ab- 
sorbed profiles. 

Although the strong absorption and blueward asymmetry 
make measurements of the widths of the UV emission lines 
difficult, we have attempted to measure the FWHMs of the 
C IV 1548.2 A and 1550.8 A lines from the STIS data. To 
do this, we started by assuming the centers of these lines are 
consistent with rest-frame emission in NGC4051. We have 
attempted to deblend the features by fitting the cubic spline 
representing the continuum in this area with two Gaussians 
centered at the fixed wavelengths of the transitions; we re- 
quired the relative strengths of the 1548.2 A and 1550.7 A lines 
to lie between 1:1 and 2:1. After obtaining the best fit, we 
blueshifted the lines by steps of 50 kms" 1 and repeated the fit- 
ting. We find the best fit to be when the lines are blueshifted by 
« 100 km s" 1 and the strengths of the lines are approximately 



equal. The FWHMs of the lines are « 1040 ±250 km s" 1 . 
These FWHMs, although rather uncertain, are consistent with 
the value of 1110± 190 km s" 1 reported by Peterson et al. 
(2000) for the FWHM of the variable broad component of H/3. 

3. DISCUSSION 

3.1. X-ray and UV Absorption and Emission 

The most significant results of this analysis are the discovery 
of X-ray emission lines and blueshifted X-ray absorption lines 
similar to those recently found from several other Seyfert Is 
(e.g., Kaastra et al. 2000; Kaspi et al. 2000), as well as the iden- 
tification of multiple intrinsic X-ray and UV absorption systems 
in the spectrum of NGC 405 1 . The velocities of the X-ray ab- 
sorption components are measured to be -2340 ± 130 km s" 1 
(the X-ray warm absorber with the highest velocity yet discov- 
ered) and -600 ±130 km s" 1 , with an average velocity for the 
emission lines (excluding Fe Ka) of 14 ± 120 km s" 1 . The low- 
velocity X-ray absorption system, taking into account its mea- 
sured width of 750^28(1 km s" 1 , is consistent with the velocity 
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range of the intrinsic UV absorption systems from -650 km s 
to -270 km s" 1 (systems 1-5). On the other hand, the high- 
velocity X-ray absorption system appears to have no kinemati- 
cally consistent counterpart in the UV (see Figure 11a), and the 
UV absorption systems with the smallest blueshifts appear to 
have no strong counterparts in the X-rays. These observations 
are physically consistent, however, because the lowest-velocity 
UV absorption systems are seen preferentially in low-ionization 
species while the high-velocity X-ray absorption system is seen 
only in very high-ionization species. Another notable point is 
that the X-ray absorption may well be subdivided into more 
than the two systems we are able to resolve with the HETGS. 
For example, Figure 1 1 shows the gain in information between 
IUE resolution and STIS resolution; there may well be equal 
complexity awaiting discovery in X-rays. This warning should 
be remembered in efforts to perform detailed modeling of the 
X-ray absorption lines, particularly when using data from the 
XMM-Newton Reflection Grating Spectrometer (RGS) which 
has lower resolution than the Chandra HETGS. 

In terms of ionization potential, the high-velocity X-ray ab- 
sorption component is observed preferentially in the higher- 
ionization (H-like) ions, while the low-velocity component is 
stronger in the lower-ionization (He-like) ions (i.e., the ioniza- 
tion potentials to create the H-like ions in Figure 3 from their 
He-like counterparts are all higher than the ionization potentials 
to create the He-like ions in this figure from their Li-like coun- 
terparts), although its presence is also statistically suggested in 
some of the H-like ions. This ionization potential approach was 
phenomenologically motivated by the observed similarities be- 
tween many of the velocity spectra for the H-like ions as a group 
and the He-like ions as a group. However, it is also relevant to 
compare the ionization parameters for which the relative pop- 



ulations of particular ions are maximized (assuming photoion- 
ization equilibrium). Kallman & McCray (1982) present curves 
showing the relative fractions of different ionization states of O, 
Ne and Si as functions of ionization parameter in the context of 
several different X-ray nebular models. Examination of these 
curves reveals that the population of He-like Si XIII is expected 
to peak at a larger ionization parameter than either of the H-like 
ions O VIII and Ne X, which seems to conflict with the ion- 
ization potential approach used to construct Figure 3. Detailed 
modeling will be required to ascertain if there is truly a conflict, 
but unfortunately such modeling will prove difficult if the ab- 
sorbing gas in NGC 405 1 is not in photoionization equilibrium 
(e.g., Krolik & Kriss (1995); McHardy et al. 1995; Nicastro 
et al. 1999). NGC 4051 typically shows large-amplitude vari- 
ability of its average flux on timescales ranging from ~ 100 s 
to a couple days (see Figure 2 and McHardy, Papadakis, & Ut- 
tley 1998), so if the recombination timescale of the warm ab- 
sorber lies in this range it is unlikely to be in photoionization 
equilibrium. All states of highly ionized oxygen, for example, 
recombine on a timescale given by f rec « 200«g 1 r 5 - 7 s, where 
n = I0 9 n 9 cm" 3 and T = 10 5 T 5 K (Shull & Van Steenberg 1982). 
Adopting a typical warm absorber temperature of 2 x 10 5 K, 
we find that the warm absorber will fall out of photoioniza- 
tion equilibrium if its density lies between ~ 10 6 cm" 3 and 
~3x 10 9 cm" 3 . 13 This fairly wide range of densities is con- 
sistent with the constraints derived using He-like lines and with 
the upper end of the range of plausible densities for the Narrow 
Line Region (see below for the He-like line constraints and dy- 
namical consistency between the Narrow Line Region lines and 
the X-ray warm absorber lines). Hence the warm absorber of 
NGC 405 1 is plausibly not in photoionization equilibrium, and 



13 If the density is <C 10 6 cm 3 , corresponding to a value of f rec much longer than the typical variability timescale, then the warm absorber might be in equilibrium 
with the average continuum flux; it would not respond to short-timescale continuum flux variations. 

Table 5 

UV Absorption Systems in the HST STIS Spectrum 









NV 


NV 


Si IV 


Si IV 


CIV 


CIV 






Velocity 


A 1238.8 


A 1242.8 


A 1393.8 


A 1402.8 


A 1548.2 


A 1550.8 


Sys. 


Velocity 


FWHM 


EW 


EW 


EW 


EW 


EW 


EW 


No. 


(kms- 1 ) 


(kms" 1 ) 


(mA) 


(mA) 


(mA) 


(mA) 


(mA) 


(mA) 


G 
1 

2 


-727 
-647 
-505 


74 
40 
165 


110+J4 
628 ±15 


125+/3 
595 ±8 


240+| 


130+H 
40+If 


301+1J 
223+|2 
750+J? 


537+^ 0a 
253+J a 
710 ±4 a 


3 


-430 


63 


123 ±7 


109 + | 






133 ±6 


100+| 


4 
5 
6 


-337 
-268 
-158 


52 
133 
45 


121+5 
425 ±7 
121 ±5 


115+1 
441+|3 
104+« 




24+11 


168+| 
537+] 0:i 
253+] a 


121+1 

418+ 6 

39+>° 


7 


-107 


64 


179 ±5 


196* 


30+J2 




710±4 a 


360+° 


8 
9 


-48 
30 


84 

23 


154+|4 




281^ 
221 


189+1° 
16+" 


710 + 4" 
152 ±3 b 


300+9 
31+' 1 



Note. — The FWHM values represent the averages from the N V lines, except for systems G and 9 which 
are from C IV and Si IV. The uncertainties on the system velocities and FWHMs are roughly 10-15 km s _1 . 
EW errors are reported at the la level. 

"Blended with one or more other components (see Figure 7 for clarification). The components that are 
blended too strongly to separate their EWs are the following: component5 (1548 A) and component G (1551 A); 
component 6 (1548 A) and component 1 (1551 A); components 7 and 8 (1548 A) and component 2 (1551 A). 

b Component 9 (1548 A) is also strongly blended with components 7 and 8 (1548 A) and component 2 
(1551 A). 
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FIG. 12. — Strongest X-ray emission lines from NGC4051. Panel (a) is from the HEG first-order spectrum, while panels (b) and (c) are from the MEG first-order 
spectrum. All features have been smoothed using a boxcar filter that is three bins in width. Note the differences in vertical scale among the panels. 



our phenomenological approach based on ionization potential 
seems appropriate for this analysis. The limited photon statis- 
tics prevent us from investigating the details of the ionization 
physics via timing analyses (e.g., Nicastro et al. 1999). 

The strongest emission features seen in the X-ray spectrum 
of NGC4051 are the O VII 22.1 A and Ne IX 13.7 A lines. 
Such features have been predicted to arise from the X-ray ab- 
sorbing gas (e.g., Krolik & Kriss 1995; Netzer 1996) and have 
now been seen from other Seyfert Is as well. The low Galac- 
tic column density to NGC 405 1 allows for a full appreciation 
of the strengths of these lines, which are shown (along with 
Fe Ka) in Figure 12. The best-fit FWHMs of these lines are 
357*220 km s" 1 and 554*320 km s -1 , respectively, which are con- 
sistent with being unresolved by the MEG. The upper limit on 
the width of the O VII line is significantly below the observed 
widths of the broad optical and ultraviolet lines. If a simple an- 
ticorrelation of line width and radial location in a virialized sys- 
tem is appropriate, this line must arise in a region farther from 
the central source than the Broad Line Region. It may arise in 
the Narrow Line Region, which for NGC 405 1 has a character- 
istic [O III] FWHM of 210-330 km s" 1 (e.g., Heckman et al. 
1981; De Robertis & Osterbrock 1984). However, we note that 
the O VII line probably arises in a higher-ionization gas compo- 
nent than the [O III] line. An alternative explanation is that this 
high-ionization line may originate in the gas that produces the 
optical coronal lines (e.g., Porquet et al. 1999). 

The Fe Ka line from NGC 4051 (see Figure 12 and Table 2) 
also deserves mention. The upper limit on the width of the 
Chandra line (FWHM) is w 2800 km s" 1 , substantially below 
the width of the broad accretion-disk line observed by ASCA 
(e.g., G96). The line is also consistent with being unresolved 
by the HEG. The line's width, along with its measured EW of 
158!$ eV mA) and energy of 6.41!°;°] keV, are con- 

sistent with an origin in neutral material surrounding the nu- 
clear region, such as the putative molecular torus of AGN uni- 



fication schemes (e.g., Krolik, Madau, & Zycki 1994). We note 
that the EW of this relatively narrow feature is a significant frac- 
tion of the EW of the line seen by G96 (« 300 eV). This point 
should be remembered in future attempts to model the broad 
Fe Ka line. 

Since the UV absorption seen from NGC 405 1 strongly ab- 
sorbs the broad UV emission-line flux, it is clear that the gas 
in which this absorption arises must be coincident with or far- 
ther from the central engine than the UV Broad Line Region. 
Our estimate of the width of the C IV lines, which is consistent 
with the width of the H/3 line, suggests that the physical sizes 
of the optical and UV Broad Line Regions are roughly compa- 
rable. This sets a rough lower limit of 3-6 light days (Peterson 
et al. 2000) on the distance from the central source to the UV- 
absorbing gas. 

3.2. Column Density and Density Estimates 

We have estimated ionic column densities using the X-ray 
absorption-line EWs and the 'curve of growth' technique de- 
scribed in § 3.4 of Spitzer (1978). In these analyses, we have 
treated the high-velocity and low-velocity X-ray absorption 
systems separately. Unfortunately, the column densities for the 
high-velocity system are quite uncertain because (1) the lines 
are not resolved so we can only constrain the velocity-spread 
parameter b to be < 1000 km s" 1 and (2) we are near or on 
the 'flat' part of the curve of growth. Furthermore, if the X-ray 
absorption lines have unresolved substructure (see Figure 11) 
this may confuse curve of growth analyses. The two X-ray ab- 
sorption lines in the high-velocity system with the best photon 
statistics are Mg XII 8.419 A and Si XIV 6.181 A. If we re- 
quire b > 100 km s" 1 , for example, the Mg XII column density 
formally allowed by our curve of growth analysis is between 



1 7 x 10 1 



and « 4 x 10 cm . We consider the upper 



end of this range to be implausibly large (e.g., the implied ab- 
sorber becomes optically thick to electron scattering), but the 
lower end is plausible. For solar abundances and the likely 



14 



COLLINGE ET AL. 



ionization correction, the corresponding lower bound on TVh is 
~2x 10 21 cm" 2 . Similarly, for Si XIV the column density for- 
mally allowed is between « 9 x 10 16 cm" 2 and « 5 x 10 20 cm" 2 . 
Assuming solar abundances, the column densities of these two 
ions could plausibly arise in a zone characterized by a single 
ionization parameter, although this is not surprising given the 
large uncertainties. We note that the O VIII column density of 
(1.4 ± 0.3) x 10 18 cm" 2 derived from the absorption edge depth 
is consistent with the plausible range of column densities in the 
high-velocity X-ray absorption line system. 

The ionic column densities of the low-velocity system are 
easier to constrain because (1) the O VII lines are resolved and 
thus we can measure b directly and (2) we are on the linear 
part of the curve of growth. The low-velocity lines most useful 
for this analysis are O VII 17.768 A, 18.628 A and 21.602 A. 
The allowed O VII column density ranges for these lines are, 
respectively, (5-20) xlO 16 cm" 2 , (3-7) xlO 16 cm" 2 , and (8- 
20)xl0 cm Although these ranges are not all consistent 
with a single value, they are not wildly different. With so- 
lar abundances and a plausible ionization correction, the im- 
plied 7V H is between <~ 1 x 10 19 cm" 2 and ~ 5 x 10 20 cm" 2 . 
We point out, however, that the O VII column density of 
(8 ± 1) x 10 17 cm" 2 indicated from the absorption-edge depth 
is larger than the upper limit derived from the curve of growth 
technique. The seeming inconsistencies among the implied col- 
umn densities may result from either the effects of iron L-shell 
absorption lines near the O VII absorption edge or the existence 
of multiple unresolved absorption components within the low- 
velocity X-ray absorption system. 

We note that the optical spectrum of NGC 405 1 contains a 
small contribution from polarized light (e.g., Grupe et al. 1998). 
If there is a corresponding scattered X-ray continuum compo- 
nent (e.g., if the flux is due to Thompson scattering), non-zero 
flux levels may be observed in the troughs of absorption lines 
that are saturated along the primary line of sight. An effect 
of this type could confuse column density estimates made us- 
ing the techniques detailed above. However, because the po- 
larization fraction is small and because no noticeable effect of 
this type is observed for the saturated UV absorption lines (see 
§ 2.2. 1), this effect should hopefully be small for the X-ray lines 
as well. 

The emission lines from the He-like ions we identified can 
serve as plasma density diagnostics (e.g., Porquet & Dubau 
2000). For both O VII and Ne IX we identified, in emission, 
only the intercombination (;) and forbidden (/) lines and not 
the resonance (r) line. The forbidden-to-intercombination line 
ratios are f/i = 4.8!^, and 5.0! 2 5 2 for O VII and Ne IX, respec- 
tively (errors are estimated by the numerical method described 
in § 1.7.3 of Lyons 1991). As the resonance emission line is not 
identified, its flux is lower than the intercombination line flux, 
and the forbidden-to-resonance line ratio (f/r) has to be > 4. 
Such high f/r implies that the plasma emitting the lines is pho- 
toionization dominated (with little collisional ionization), and 
the upper limit on its temperature is 10 6 K. The upper limit on 
the plasma density, using the f/i ratios and Figure 8 of Porquet 
& Dubau (2000), is then 4 x 10 10 cm" 3 . 

The S/N of the X-ray spectrum of NGC 4051 is not fa- 
vorable for detailed photoionization modeling such as Kaspi 
et al. (2001) have performed for NGC 3783. However, as Fig- 
ure 13 shows, the correlation between the EWs of spectral fea- 
tures from NGC 405 1 and from NGC 3783 is striking (although 
NGC 3783 has stronger X-ray absorption overall). This sug- 



gests that the basic physical nature of the absorbing gas is sim- 
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FIG. 13. — Comparison of the EWs of features detected in the Chandra 
spectra of both NGC 4051 and NGC 3783 (Kaspi et al. 2000). Dots repre- 
sent absorption features, and triangles represent emission features. EWs of 
absorption lines from NGC 4051 that are detected in both the high- velocity 
and low-velocity systems represent the sums of the measured values for the 
two systems. A line with slope equal to unity is plotted to guide the eye. Al- 
though the features seen from NGC 405 1 are generally weaker than those seen 
from NGC 3783, the noticeable correlation spanning a factor of 100 in EW 
indicates a basic similarity in the warm absorbers of the two systems. 

ilar in these two Seyferts and that the results of the detailed 
modeling of the absorber in NGC 3783 should be applicable, at 
least to first order, for NGC 4051 as well. This correlation is 
notable given (1) the different 0.5-10 keV ionizing continuum 
shapes of the two objects (while their average hard power-law 
slopes are similar, NGC 3783 lacks the strong soft X-ray excess 
seen from the NLS 1 NGC 405 1) and (2) the fact that NGC 405 1 
is w 100 times less luminous in X-rays than NGC 3783. How- 
ever, it is also important to note that the UV absorption in 
NGC 4051 is much stronger than in NGC 3783 (while the op- 
posite is true for the X-ray absorption). Crenshaw et al. (1999) 
report EWs of 0.28 ±0.04 A and 0.19 ±0.03 A for the N V 
1238.8 A and 1242.8 A lines from NGC 3783, while the total 
values we measure for NGC 405 1 are roughly seven and nine 
times larger, respectively. The combination of weaker X-ray 
absorption and stronger UV absorption seems to indicate that 
the absorbing gas in NGC 405 1 is characterized by a smaller 
ionization parameter than the gas in NGC 3783. This may result 
from the different extreme-ultraviolet to soft X-ray continuum 
shapes of the two objects (since this is the ionizing continuum 
relevant for the UV lines) or from the difference in luminosity 
between the two objects; observations of additional Seyferts are 
required to search for connections of these types. 

3.3. Variability and Soft X-ray Excess 

The low state that NGC 4051 entered for the last w 15 ks of 
our Chandra observation is comparable to the lowest states in 
which it has been previously observed, although it is not the 
lowest. Guainazzi et al. (1998) reported a state in which the av- 
erage 2-10 keV flux was 1.3 x 10~ 12 ergcm~ 2 s" 1 , smaller by a 
factor of about four than that which we observe. The flux during 
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our low state is comparable to that seen during the 1990 Broad- 
Band X-ray Telescope (BBXRT) low-state observation (Weaver 
1993). Previous analyses have found that the 2-10 keV spec- 
trum of NGC 405 1 becomes increasingly hard with decreasing 
flux (e.g., Matsuoka et al. 1990; G96; Uttley et al. 1998). Our 
data show the expected qualitative behavior (see §2.1.3), and 
quantitatively our hard photon index of T = 0.83+y [9 for the low 
state is consistent with the T = 0.8 ±0.3 observed by BBXRT (as 
well as the T = 0.78^;^ observed by Guainazzi et al. 1998). 
However, both of these values lie somewhat below a linear ex- 
trapolation of the photon index versus X-ray flux correlation 
established for NGC 4051 at higher fluxes; either this correla- 
tion becomes nonlinear at low fluxes or we are witnessing the 
effects of its long-term secular evolution. The flat photon index 
is probably due to a large contribution from a reflection com- 
ponent that lags the primary continuum variations, but unfortu- 
nately our ps 15 ks of low-state Chandra data lack the photon 
statistics to examine this matter in detail. 

Our high-resolution data on the soft X-ray excess reveal that 
(1) it is continuous rather than composed of a forest of narrow 
emission lines, (2) it shows significant spectral curvature, and 
(3) it is rapidly variable. A simple blackbody model does not 
entirely fit the shape of this excess, and comparison with more 
detailed physical models is required. The continuous nature of 
the soft X-ray excess is as expected given its rapid variability 
(Elvis et al. 1991), and it parallels the recent result of Turner 
et al. (2001) regarding the NLS1 Ton SI 80. 

4. CONCLUSIONS 

We have performed the first X-ray grating spectroscopy and 
the first detailed UV absorption study of the famous Narrow- 
Line Seyfert 1 galaxy NGC 405 1 . Our high-resolution X-ray 
and UV spectroscopic analyses reveal a number of significant 
results: 

• We detect a number of O, Ne, Mg and Si X-ray ab- 
sorption lines belonging to two distinct absorption sys- 
tems at blueshifted velocities of -2340 ± 130 km s" 1 
and -600 ± 130 km s" 1 . The former, which has the 
highest velocity of any known X-ray warm absorber 
in a Seyfert, is observed preferentially in the higher- 
ionization (H-like) ions, while the latter is stronger in 
lower-ionization (He-like) ions (see § 2.1.2, § 3.1, Ta- 
bles 1 and 2, and Figures 1 and 3). 

• The STIS spectrum reveals rich intrinsic UV absorption 
consisting of multiple (ps 10) systems, some of which 
are consistent in velocity with the low-velocity X-ray 
absorption system (but not the high-velocity system) 
and some of which are detected in both high-ionization 
(C IV, N V and Si IV) and low-ionization (C II, Si III and 
Si II) ions. Our analysis of archival IUE data reveals 
that UV absorption has been present at some level for 



more than 20 yr (see § 2.2.2, § 2.2.3, §3.1, Table 5, and 
Figures 6, 7, 8 and 10). 

• The low-velocity X-ray absorption system is resolved 
by the MEG. We find it to have a FWHM of 
750^280 km s" 1 . Comparison with the UV absorption in 
the same velocity range suggests that the X-ray absorp- 
tion systems may have complex, unresolved substruc- 
ture (see § 2.1.2, § 3.1, and Figures 4 and 11). 

• The Chandra spectrum contains several unresolved X- 
ray emission lines, the strongest of which are from O VII 
and Ne IX. These lines have velocities consistent with 
rest-frame emission in NGC 4051 (see § 2.1.2, § 3.1, 
Table 2, and Figures 1 and 12). 

• We estimate that the X-ray warm absorber has a col- 
umn density of Nh ~ 10 21 cm -2 in the high- velocity 
system and V H ~ 10 20 cm" 2 in the low- velocity system 
(although we note that these values may be subject to 
systematic uncertainties). Plasma diagnostic techniques 
reveal that the gas is primarily photoionized (rather than 
collisionally ionized) and set upper limits of ps 10 6 K 
and ps 4 x 10 10 cm" 3 on its temperature and density, re- 
spectively (see § 3.2). 

• We find that the soft X-ray excess is rapidly variable, 
that it is continuous rather than made up of many nar- 
row emission lines, and that it has significant spectral 
curvature (see § 2.1.3, § 3.3, Table 4, and Figure 5). 

• During the last ps 15 ks of our Chandra observation, 
NGC 4051 entered a low-flux state comparable to the 
lowest states in which it has been previously observed. 
In the low state the spectrum became substantially 
harder, as in previously documented cases (see § 2.1.3, 
§ 3.3, and Figures 2 and 5). 
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